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This paper presents updated limits on "disappearance 
decays" of B° mesons [l[ , where the B° decay contains no 
observable final state particles, or such "invisible" decay 
products plus a single photon. We define invisible decay 
products here to be electrically neutral particles that do 
not generate a signal in the electromagnetic calorimeter. 
These results represent an improvement over the previous 
limits on these decays, which were based on 19% of the 
present data sample [2j. 

The rate for invisible B decays is negligibly small 
within the Standard Model (SM) of particle physics, but 
can be larger in several models of new physics. The SM 
decay B — > w, which would give such an invisible ex- 
perimental signature, is strongly helicity-suppressed by a 
factor of order (m^/m^o) 2 [3[ and the resulting branch- 
ing fraction is necessarily well below the range of present 
experimental observability. The SM expectation for the 



B° — > vv^j branching fraction is predicted to be of order 
10~ 9 , with very little uncertainty from hadronic inter- 
actions [4] . An experimental observation of an invisible 
(+ 7) decay of a B° with current experimental sensitivity 
would thus be a clear sign of physics beyond the SM. 

A phenomenological model motivated by the observa- 
tion of an anomalous number of dimuon events by the 
NuTeV experiment 5] allows for an invisible B° decay 
to a Pxl final state, where x? is a neutralino, with a 
branching fraction in the 10~ 7 to 10~ 6 range [6|. Also, 
models with large extra dimensions, which would provide 
a possible solution to the hierarchy problem, can have the 
effect of producing significant, although small, rates for 
invisible B° decays [7H9J. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II e+e~ collider at SLAC. 
The data sample corresponds to a luminosity of 424 fb _1 
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accumulated at the T(4S) resonance and contains (471 ± 
3) x 10 6 BB pair events. For background studies we also 
used 45 fb _1 collected at a center-of-mass (CM) energy 
about 40 MeV below BB threshold (off-peak). 

A detailed description of the BABAR detector is pre- 
sented in Ref. [ltj. Charged particle momenta are mea- 
sured in a tracking system consisting of a five-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
hexagonal-cell wire drift chamber (DCH) . The SVT and 
DCH operate within a 1.5 T solenoidal field, and have a 
combined solid angle coverage in the CM frame of 90.5%. 
Photons and long-lived neutral hadrons are detected and 
their energies are measured in a CsI(Tl) electromagnetic 
calorimeter (EMC), which has a solid angle coverage in 
the CM frame of 90.9%. Muons are identified in the in- 
strumented flux return (IFR). A detector of internally 
reflected Cherenkov light (DIRC) is used for identifica- 
tion of charged kaons and pions. A GEANT4 [ll[ based 
Monte Carlo (MC) simulation of the BABAR detector re- 
sponse is used to optimize the signal selection criteria 
and evaluate the signal detection efficiency. 

The detection of invisible B decays uses the fact that 
B mesons are created in pairs, due to flavor conserva- 
tion in e + e~ interactions. If one B is reconstructed in an 
event, one can thus infer that another B has been pro- 
duced. We reconstruct events in which a B° decays to 
]j(*)~l+ v (referred to as the "tag side"), then look for 
consistency with an invisible decay or a decay to a single 
photon of the other neutral B (referred to as the "signal 
side" ) . The choice of reconstructing semileptonic B° de- 
cays on the tag side, with respect to fully-reconstructed 
B° final states, is motivated by a higher reconstruction 
efficiency. A disadvantage is the presence of the invisible 
neutrino, which prevents the exploitation of kinematic 
variables such as the reconstructed B° mass. However, 
the background contamination is mitigated by the pres- 
ence of a high momentum lepton. 

We reconstruct D*~ mesons in the final states D°tt~ 
or D~tt°, with D° decays to K + ir~, K + ir~ir Q , or 
K + ir~ir + Tr~ , and D~ decays to K + ir~ir~ or K®ir~ . 
We identify K + candidates using Cherenkov-light in- 
formation from the DIRC and energy-loss information 
(dE/dx) from the DCH and SVT. The K° mesons are 
reconstructed in the decay mode K® —¥ ir + ir~, where 
the 7r + 7r~ invariant mass lies in a ±25 MeV/c 2 win- 
dow around the nominal Particle Data Group (PDG) K® 
mass [12| . The ir° candidates are composed of pairs of 
photons observed in the EMC. Each photon must have 
a reconstructed energy above 30 MeV in the laboratory 
frame, and the sum of their energies must be greater than 
200 MeV. The ir° candidates must have an invariant 
mass between 115 and 150 MeV/c 2 . A mass-constrained 
fit is imposed on ir° candidates in order to improve the 
resolution on the reconstructed invariant mass of the par- 
ent D meson. 

Kaon and pion candidates are then combined to recon- 



struct mesons. These are required to have an invari- 
ant mass within 60 MeV/c 2 of their nominal PDG mass, 
except for D decays with a ir° daughter, which must 
be within 100 MeV/c 2 of the nominal D° mass. Mass- 
constrained fits are applied to D° and D~ candidates in 
order to improve the measurement of the momentum of 
each D. The difference in reconstructed mass between 
D*~ decay candidates and their D daughters must be 
less than 175 MeV/c 2 and greater than 137 MeV/c 2 . All 
candidates must have a total momentum between 
0.5 and 2.5 GeV/c in the CM frame. 

Tracks selected as lepton candidates must pass either 
electron or muon selection criteria. We identify electron 
candidates using energy and cluster shape information 
from the EMC, and Cherenkov angle information from 
the DIRC. Muon candidates are identified using informa- 
tion from the IFR and EMC. Both electrons and muons 
must also have a momentum of at least 0.8 GeV/c in the 
laboratory frame, and a minimum of 20 DCH measure- 
ments. 

To further select B° — > i + v candidates, we re- 

quire a Z}(*) _ candidate and a lepton candidate to be 
consistent with production at a common point in space. 
The decay vertex is reconstructed from a kinematic fit 
to all the candidate daughters, and a minimum \ 2 ver ~ 
tex probability of 0.001 is required. We then calculate 
the cosine of the angle between the D^~£ + and the hy- 
pothesized B° candidate in the CM frame, under the 
assumption that the only particle missing is a neutrino: 



cos 9 



2 EbE D (,)- 1+ — m B 



B,D(*)-i+ 



2 \pb\\p D (*)-i+\ 



(1) 



The energy in the CM frame E D (,)-^+ and mass to £><*)-£+ 
of the Z)(*)~ £ + combination are determined from recon- 
structed momentum information, and tub is the nominal 
B° mass [l^ . The B° momentum \pb\ and energy Eb in 
the CM frame are determined from beam parameters. If 
our assumption that there is only one missing particle, a 
neutrino, in the B° decays is incorrect, cos 6g DW—&+ can 
fall outside the region [—1,1]. We require the £>W-£+ 
combination to satisfy —5.5 < cos B d(*)-i+ < 1-5- The 
selected region allows for non-physical cos9 B r>(.*)-e+ val- 
ues, accounting for detector energy and momentum res- 
olution. Moreover the asymmetric cut admits higher D* 
mass states where additional decay products are lost. In 
the rest of the analysis such products are not associ- 
ated with the tag side decay chain but are considered 
as extra particles in the event. When more than one 
£>° — > Z)W _ i + v candidate is reconstructed in an event, 
the one with the highest vertex probability is taken. 

We consider events with no charged tracks besides 
those from the B° — !• D^~£ + v candidate. In order 
to reject background events where one charged or neu- 
tral particle is lost along the beam pipe, the cosine of 
the polar angle of the missing momentum in the CM 
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FIG. 1: Distributions of the NN output for simulated B° — > 
invisible events with a D meson on the tag side. The black 
solid line is the signal while the red dashed line is the back- 
ground. The solid gray vertical line defines the NN output 
signal region. 



frame (cos 0^ ss ) is required to lie in the [—0.9, 0.9] range. 
The missing 4-momentum due to unreconstructed parti- 
cles is defined as the difference between the T(AS) and 
the reconstructed tag side 4-momentum. fn the B° — > 
invisible+7 channel the signal-side photon 4-momentum 
is also subtracted from the T(AS) one. 

For the B° — > invisible decay, in events where the D 
meson on the tag side decays into K^n+ir^, two addi- 
tional selection criteria are applied. The first concerns 
the sum of the cosine of the angles between the kaon and 
two pions, cosOktt! + cos0ktt 2 > —0.8, while the second 
concerns the sum of the cosine of the angles between the 
lcpton and the pions, cos6^ 1 + cos#£ 7r2 < 0.8. The main 
effect of this selection is the reduction of the background 
from misreconstructed e + e~ — » t + t~ events. 

To reconstruct B° — > invisible+7 events, one remain- 
ing photon candidate with energy greater than 1.2 GeV 
in the CM frame is also required. If the detected photon 
has an energy smaller than 1.2 GeV in the CM frame, the 
event falls in the B° — > invisible category and the neutral 
candidate is considered as an extra photon in the event. 
The choice of this cut generates a cross-feed between the 
two channels; MC simulation studies show that this has 
a negligible effect on the final result. 

A neural network (NN) is used to provide further 
discrimination between signal and background events. 
We use the TMVA software package [13| and its mul- 
tilayer perceptron implementation of a NN. The NN is 
trained using samples that represent the signal and back- 
ground components. For the signal, a MC simulated 
sample is used, in which a generic semileptonic B decay 
is generated and reconstructed. Weighted off-peak data 
(composed of e + e~ — > ce,uu,ss, and t + t~ events, de- 
noted as continuum background) and MC simulated BB 



events are used to describe the background contamina- 
tion. Off-peak data are used to model continuum back- 
ground, as the MC was found to incorrectly reproduce 
the cross-section of two-photon fusion events, such as 
e + e~ — > e + e~77 — > e + e~qq/r + T~. These events typ- 
ically have decay products directed along the beam lines, 
and thus outside the detector acceptance. 

The variables used as input for the NN, common to the 
B° — > invisible and B° — > invisible+7 analyses, are: 1) 
cos 9 B £>(*)-£+ ; 2) the cosine of the angle in the CM frame 
between the thrust axis (the axis along which the total 
longitudinal momentum of the event is maximized) and 
the flW"^ pair momentum direction; and 3) the lep- 
ton momentum in the CM frame. In the B° — > invisible 
analysis, we additionally use 1') M^f ss (defined as the 
invariant mass of the event after the D^~£ + pair is sub- 
tracted); 2') the B meson vertex fit probability; 3') the 
ratio between the first and the zeroth-order L momenta 
in the CM frame: 



Li = ^^pcos 1 6, 



(2) 



where the sum is over extra tracks and neutrals and is 
computed with respect to the thrust axis; 4') the trans- 
verse momentum of the D~l + pair in the CM frame; 5') 
the minimum invariant mass of any two charged tracks 
in the event; and 6') the minimum invariant mass of any 
three charged tracks in the event. Variables 4') - 6') enter 
the NN only in the case of a reconstructed B° — > D~ ' l + v 
decay on the tag side. In the B° — > invisible+7 analy- 
sis, wc additionally use 1") the energy of the photon on 
the signal side evaluated in the laboratory frame and 2") 
Mmiss (f° r B° — > D~l + v reconstructed events only). 

The selection on the output of the NN is optimized 
by minimizing the expected upper limit on the branch- 
ing fraction, defined by a Bayesian approach as detailed 
later in this paper, under the hypothesis of observing zero 
signal events. In Fig. [TJ the output of the NN for simu- 
lated B° — > invisible with a D meson on the tag side and 
the corresponding signal region are shown. 

After the NN selection, the D meson invariant mass 
(tud) and the difference between the reconstructed D* 
invariant mass and the PDG D° mass (Am) are used to 
define a signal region (SR) and a side band region (SB) 
for the D tag and D* tag samples, respectively. The SR 
is defined as a ±15 MeV/c 2 window around the PDG 
value for run for the B° — > D~l + v sample, and as 0.139 
< Am < 0.148 MeV/c 2 for the B° -> D*~l + v sample. 
The excluded regions are used as the SB region. 

The total energy in the EMC computed in the CM 
frame and not associated with neutral particles or 
charged tracks used in the D^*^£ + reconstruction is de- 
noted as -Eoxtra- For B Q — > invisible+7, the energy of the 
highest-energy photon remaining in the event (the signal 
photon candidate) is also removed from the -Eextra com- 
putation. The -Ecxtra signal region is defined by impos- 
ing an upper bound at 1.2 GeV. In both B° — > invisible 
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and B° — > invisible+7 samples, this variable is strongly 
peaked near zero for signal, whereas for the background 
the distribution increases uniformly in the chosen signal 
region. Background events can, however, populate the 
low -Ecxtra region, when charged or neutral particles from 
the event are either outside the fiducial volume of the 
detector or are unreconstructed due to detector ineffi- 
ciencies. Contributions from misreconstructed tt° decays 
usually populate the high E ex tra region. 

Using detailed Monte Carlo simulations of B° — > in- 
visible and B° — > invisible+7 events, we determine our 
signal efficiency to be (17.8 ±0.2) x 1CT 4 for B° -> invis- 
ible and (16.0 ±0.2) x 10~ 4 for B° — > invisible+7, where 
the uncertainties are statistical. 

We construct probability density functions (PDFs) for 
the -Eextra distribution for signal (V s ig) and background 
('Pbkg) using detailed MC simulation for signal and data 
from the mn and Am sidebands for background. The 
two PDFs are combined into an extended maximum like- 
lihood function C, defined as a function of the free param- 
eters 7V s ig and iVbkg, the number of signal and background 
events, respectively: 



C(N sig ,N hkg 



[(1 - z sig )iV sig + (1 - z hkg )N hkg ] Nl 



x e -[(l-2=i g )^i g + (i-Zbk g )AW] 
T^sig (-^cxtra,z |Psi 



n 



7 , bkg(£ ; oxtra,i IPbkg) 



(1 - Z sig )N si 

m 

(1 - Zbkg)AW 



(z sig N sig + Zbkg^Vbkg) 



N 



x e 



N ! 

— (•Zsig-N s ig + ZbkgA r bkg) 



(3) 



The photon reconstruction has a detection lower energy 
bound of 30 MeV, and as a consequence, the -Eextra dis- 
tribution is not continuous. To account for this effect, 
the likelihood in Eq. [3]is composed of two distinct parts, 
one for -Eextra > 30 MeV and one for E cx tra = MeV. In 
the likelihood function, z s i g and Zbkg are the fractions of 
events with E ex tra = MeV for signal and background, 
respectively, and p s - lg and pbkg are the vectors of parame- 
ters describing the signal and background PDFs, a kernel- 
based PDF |14| and a second-order polynomial, respec- 
tively. The fixed parameters No, N±, and E ex tra,i are, 
respectively, the number of events with E ex tra = MeV, 
the number of events with -E ex tra > 30 MeV, and the 
value of Eextra for the i th event. 

The negative log-likelihood is then minimized with re- 
spect to -V s jg and Vbkg in the data sample. The resulting 
fitted values for N s i g and A^kg are given in Table |TJ Fig- 
ure shows the E ex tra distributions for B° — > invisible 
and B° — > invisible+7 with the fit superimposed. 

The fitted signal yields are used to determine the decay 
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FIG. 2: Results of the maximum likelihood fit of -Eextra for 
B° — > invisible (left) and B° — s> invisible+7 (right). 



branching fractions (£>), which are defined as: 



B = 



N 



Slg 



e x N BS 



(4) 



where e is the total signal efficiency, corrected for data- 
MC discrepancies (described below), and N B § is the 
number of produced BB pairs. 

The systematic uncertainty on the signal efficiency is 
dominated by data-MC discrepancies in the distribution 
of the variables used as input to the NN. This results in 
relative uncertainties of 6.1% and 8.2% for B° — > invisible 
and B° — > invisible+7, respectively. This uncertainty is 
evaluated using the hypothesis that the data-MC agree- 
ment could reduce the discriminating power of each input 
variable. In order to make the signal distributions more 
background-like, in the signal sample we apply a Gaus- 
sian smearing to each of the NN input variables, where 
the smearing parameters are evaluated by comparing the 
difference in the root- mean-square of the signal and back- 
ground shapes. The NN output selection is then applied 
to this new sample and the difference between the nomi- 
nal signal efficiency and this new efficiency is used as the 
systematic uncertainty. 

Another important contribution is due to the estima- 
tion of the efficiency on the tag side reconstruction (3.5% 
for both channels). For this purpose, data and MC sam- 
ples in which a B° and a B° are both reconstructed as 
decays to D^£u in the same event ("double tag" events) 



TABLE I: Fitted yields of signal and background events in 
data. The uncertainties are statistical. 



Mode 




-Vbkg 


B° — > invisible 


-22 ± 9 


334 ± 21 


B° —>■ invisible+7 


-3.1 + 5.2 


113 + 12 
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are used. The square root of the ratio between the num- 
ber of the selected double tag events in data and in MC 
simulation is 0.928 (0.824) for events with B a -> D^Hv 
on the tag side; these ratios are used to correct the effi- 
ciency. The propagation of the statistical errors on the 
correction factors is used as a systematic uncertainty on 
the signal efficiency. 

Other contributions to the systematic uncertainty on 
the signal efficiency come from the choice of the preselec- 
tion criteria and from the SR definition of TO£>(Ato). The 
first effect is evaluated by applying a Gaussian smearing 
to the variables involved (cos 0^ niss , cos Qktx x + cos 0Kir 2 
and cos6g 7Tl + cos#&r 2 ). The variation on the signal effi- 
ciency is then used as a systematic uncertainty. As was 
done for the NN, this uncertainty is evaluated using the 
hypothesis that the discrimination power of each variable 
is reduced. The second effect is evaluated by changing 
each of the bounds of the SR definition by a value 6 
(3 MeV for mo and 1.5 MeV for Am), which is half of 
the me/Am resolution as evaluated in data. The rel- 
ative maximum variation in efficiency is then used as a 
systematic uncertainty. 

An additional source of systematic uncertainty is de- 
termined for the B° — > invisible+7 decay in order to 
account for detector inefficiency in the single photon re- 
construction. This is evaluated by comparing the data 
and MC ir° reconstruction efficiency in r — > p(ir ir°)v 
decays, where the total number of produced tt° in the se- 
lected sample is determined from the branching fraction 
of the specific r decay [13]. Then the ratio between the 
two efficiencies, combined with the error on the r decay 
branching ratio, is used to extract a systematic error for 
the single photon reconstruction efficiency. 

The total systematic uncertainty on the signal selection 
efficiency is 7.7% for B° —5- invisible decay and 9.5% for 
B° — »■ invisible+7 decay. 

The systematic uncertainty on the number of signal 
events is dominated by the parametrization of the back- 
ground -Eoxtra distribution. A maximum likelihood fit of 
-E-cxtra with the background parameters varied according 
to their statistical error and correlations is performed. 
For each parameter the difference in the fitted signal 
yield with respect to the nominal value is used as a sys- 
tematic uncertainty. Other contributions to the signal 
yield systematic uncertainty come from the signal shape 
parametrization and from the use of the data SB for the 
determination of the background shape. The first is eval- 
uated as the difference between the fitted yield with the 
polynomial shape and an alternative exponential shape. 
The latter, computed as the difference in the £" C xtra shape 
between the SR and SB, is parametrized with a first- 
order polynomial using the charge-conservation violating 
B + — > invisible(+7) control sample discussed below. This 
parametrization is used to weight the background shape, 
and the difference in the fitted yield is used as a sys- 
tematic uncertainty. Another contribution for the B° — > 



TABLE II: Summary of the systematic uncertainties. 



Source B° 


— > invisible B - 


■¥ invisible+7 


Normalization Errors 


B-counting 


0.6% 


0.6% 


Efficiency Errors 


Tagging Efficiency 


3.5% 


3.5% 


mo (Am) Selection 


1% 


1.3% 


Preselection 


3% 


2.4% 


Neural Network 


6.1% 


8.2% 


Single Photon 




1.8% 


TOTAL 


7.7% 


9.5% 


Yield Errors (events) 


Background Param. 


15.8 


6.5 


Signal Param. 


2.0 


1.2 


Fit Technique 




1.0 


-Eoxtra Shape 


0.1 


1.8 


TOTAL 


15.9 


6.9 



invisible+7 decay is due to a small bias observed in MC 
studies of the yield extraction. The total systematic un- 
certainties on the signal yield are 16 and 7 events for 
B° — > invisible and B° — > invisible+7, respectively. 

For the systematic contribution due to the uncertainty 
on the estimation of the total number of BB events in 
the data sample, the procedure adopted is described in 
Ref. [Hj] and the resulting uncertainty is 0.6%. The sys- 
tematic uncertainties are summarized in Table ITT1 

A Bayesian approach is used to set 90% confidence level 
(CL) upper limits on the branching fractions for B° — > 
invisible and B° — > invisible+7. Flat prior probabilities 
are assumed for positive values of both branching frac- 
tions. Gaussian likelihoods are adopted for signal yields. 
The Gaussian widths are fixed to the sum in quadra- 
ture of the statistical and systematic yield errors. We 
extract a posterior PDF using Bayes' theorem, including 
in the calculation the effect of systematic uncertainties 
associated with the efficiencies and the normalizations, 
modeled by Gaussian PDFs. Given the observed yields 
in Table HI the 90% confidence level upper limits are cal- 
culated, after the marginalization of the posterior PDF, 
by 

pUL I poo 

J V{B)dB J V(B)dB = 0.9. (5) 

The resulting upper limits on the branching fractions are 

B(B° -> invisible) < 2.4 x 10~ 5 
B(B° -> invisible + 7) < 1.7 x 10~ 5 

at 90% CL. In order to cross-check the results of the anal- 
ysis, we also search for the charge-conservation violating 
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modes B + — > invisible and B + — > invisible+7. We check 
that their resulting signal is consistent with zero. For 
these modes, we reconstruct B^- — > D°£vX°, where X° 
can be a photon, ir°, or nothing. The D° is reconstructed 
in the same three decay modes as in B° ->. £>(*)-£+ 1/ ; and 
similar criteria are enforced for the reconstructed charged 
B as for the neutral B modes. The resulting fitted val- 
ues of iV s ig are —4.3 ±3.8 (stat.) for B + — > invisible and 
-7.9 ± 8.3 (stat.) for B + invisible+7, which are both 
consistent with zero within 1.1 standard deviations. 

In summary, we obtain improved limits on branching 
fractions for B° decays to an invisible final state and for 
B° decays to invisible+7. The upper limits at 90% con- 
fidence level are 2.4 x 10~ 5 and 1.7 x 1CU D for the B a — > 
invisible and B° — > invisible+7 branching fractions, re- 
spectively. The latter limit assumes a photon momentum 
distribution predicted by the constituent quark model for 
B° — > vv^f decay whereas the B° — > invisible limit 
is not decay- model dependent. These limits supersede 
our earlier results , which used a small fraction of our 
present dataset. 
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